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We study nonlinear dynamics of the kicked particle whose motion is confined by square billiard. 
The kick source is considered as localized at the center of square with central symmetric spatial 
distribution. It is found that ensemble averaged energy of the particle diffusively grows as a function 
of time. This growth is much more extensive than that of kicked rotor energy. It is shown that 
momentum transfer distribution in kicked billiard is considerably different than that for kicked free 
particle. Time-dependence of the ensemble averaged energy for different localizations of the kick 
source is also explored. It is found that changing of localization doesn't lead to crucial changes in 
the time-dependence of the energy. Also, escape and transport of particles are studied by considering 
kicked open billiard with one and three holes, respectively. It is found that for the open billiard with 
one hole the number of (non-interacting) billiard particles decreases according to exponential law. 



I. INTRODUCTION 

Periodically driven dynamical systems play one of 
the central role in classical and quantum chaos theory 
[l[-0. An important feature of periodically driven dy- 
namical systems is the chaotization of the motion under 
certain conditions (resonances, exceeding critical value of 
the external field strength etc.). This chaotization leads 
to the exponential divergence of neighboring trajectories 
in the phase space and diffusive growth of the energy of 
the given system 

Comprehensive theoretical [H-jsl and experimental Q 
study of the simplest periodically driven systems kicked 
rotor shows that for higher enough values of the kicking 
force the average energy of the system grows linearly in 
time. Depending on the kicking strength dynamics of the 
system can be mixed or chaotic [H-ja]. 

In this paper we study particle motion in closed and 
open square billiards in the presence of external periodic 
perturbation with the aim to explore the effect of con- 
finement on periodically driven dynamics. The motion 
of a particle in confined geometries is a paradigm for the 
study of nonlinear dynamics and chaos in theoretical [5j- 
[l3l | and experimental contexts 14 1- 16 1 in case of time- 
dependent systems. Up to now much progress is done 
in the study of billiards with static boundaries. Also, 
classical dynamics of the billiards with time- dependent 
boundaries have been extensively studies in the context of 
Fermi acceleration l r il-\2~h. In particular, possibility of 
Fermi acceleration in Sinai fUj, stadium and ellip- 



geometry of its boundaries, while kicked billiard has addi- 
tional factors which allows to manipulate by the dynam- 
ics. Those are perturbation parameters such as coupling 
constant and frequency of the kick. 

We explore the kicked billiard particle dynamics by cal- 
culating time dependence of the energy both for a single 
trajectory and ensemble of the trajectories. In addition, 
we treat momentum transfer distribution for this system 
and compare it with that of kicked rotor. 

Our study shows that depending on the type of the kick 
potential, localization and strength of the perturbation 
force the dynamics can be different. In particular, the 
motion of the particle can be localized (trapped) in the 
kicking area. The 'life time' of such trap depends on 
the perturbation strength (coupling constant) and initial 
energy of the particle. 

It is found that the average energy of the kicked parti- 
cle grows diffusively as a function of time, as does the 
kicked rotor average energy. However, this growth is 
more rapid than that of kicked rotor. 

Also, we study particle transport in kicked billiard by 
considering open billiards with one and three holes. We 
explore time-dependence of the escape rate and transmis- 
sion and reflection coefficients. We note that classical dy- 
namics of unperturbed open billiards have been studied 
for integrable and chaotic geometries by many authors 
(see, e.g. 



Refs. 



24|-[3 



In particular, it was found 



tical \2u , \2i] billiards has been shown. 

Dynamics of unperturbed billiard is governed by the 



that the number of (non-interacting) particles in non- 
integrable open billiard decreases exponentially, while in 
case of regular billiard it decreases according to power 
law 24j-[26j. In this paper we extent these studies for 



the case when non-interacting particles in an open bil- 
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liard subjected into the influence of delta-kicking force. 

The motivation for the study of periodically driven bil- 
liards is caused by several reasons. In many systems 
with confinement (e.g. quantum dots, graphene, MIT 
bag model) in real situations a confined system is sub- 
jected to the action of external time-dependent fields. 
For example, hadrons in quark-gluon plasma can be con- 
sidered as driven confined system. In quantum dots or 
in graphene, external time-periodic perturbation can be 
used for manipulating by particle transport in these sys- 
tems. 

This paper is organized as follows. In the next section 
we will give formulation of the problem for the kicked 
square billiard and its detailed solution. Section III ex- 
tends the results of section II for the case of kicked open 
billiard with one and three holes. Last section provides 
some concluding remarks. 



where p and p' are the momenta just before and right 
after the kick respectively. In Fig. [3J the energy E(t) 








FIG. 1. (Color online) Three-dimensional plot of the kicking 
potential for a = 1.0. 



II. KICKED BILLIARD VS KICKED ROTOR 

In this work, we examine a system consisting of a parti- 
cle moving inside a two-dimensional square billiard with 
an additional kicking source located at the center of the 
billiard. The kicking potential (see Fig. Q} is given by 



V(x,y,t) 
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where a,T, a and R are the coupling constant, the kick- 
ing period, the side length of the square and the radius of 
kicking area, respectively. Furthermore we fix R = a/4. 
The Hamiltonian of the system can be written as 



H = H + V(p,t), 



(2) 



with Ho being the Hamiltonian of the particle moving in 
the square billiard without any kicking. Since particles 
move balistically in between collision with the billiard 
boundary as well as in between kicks, a discrete mapping 
is used to solve the corresponding equations of motion. 
At collisions with the billiard boundary the particles un- 
dergo elastic reflections, whereas if the particles are at 
t = nT, n = 1,2,3,... inside the kicking area, their mo- 
menta changes according to 




p' = p- VV(x,y,t), 



(3) 



FIG. 2. (Color online) Time-dependence of the energy for a 
typical particle in the kicked billiard (a = 0.2, T = 0.01). 
Inset: rapid oscillations of the energy, the particle is trapped 
for a certain time inside the kicking area. 

as a function of time for a typical trajectory is shown. 
The curve can be decomposed in characteristic segments, 
each corresponding to a representative dynamics of the 
particle: Firstly, parts with rapid oscillations of the en- 
ergy can be seen, see the inset of Fig. [2] During such 
phases, the particle is trapped for certain time inside the 
kicking area and experiences successive kicks, leading to 
the typical oscillations of the energy shown in the inset 
of Fig. [21 Secondly, single vertical lines in the curve of 
Fig. [2] correspond to single kicks, the energy of the parti- 
cle is either increased or decreased, depending on where 
the kick happens. Finally, parts of the curve where the 
energy remains constant. During such times, the particle 
moves outside the kicking area (quasiperiodic orbits) or 
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FIG. 3. (Color online) Distribution of constant energy time 
intervals. Circles are the results of direct numerical calcula- 
tions, solid line is the linear regression 



crosses the kicking area without getting a kick, the latter 
means it is in a way not synchronized with the kicking 
period. 

To understand the kicked billiard dynamics more 
deeply we need to explore distribution, N(t) of the time 
intervals during which the energy remains constant, i.e. 
constant-energy-time intervals. Fig. 3 presents the plot 
of such distribution in double logarithmic scale obtained 
using direct numerical computation. In addition, this 
figure compares also N(t) with the curve N(t) obtained 
from the linear regression. Such power law (with the ex- 
ponent equal to —3) behavior of N(t) can be explained 
as follows: appearing constant energy intervals is caused 
by two factors. One of them corresponds to the situa- 
tion when billiard particle moves along the quasiperiodic 
orbits which initially doesn't cross the kicking area. For 
this case one can obtain the estimate N(t) ~ t -1 . How- 
ever, we found that the probability for appearing such 
orbits in our system is quite small. Therefore such regime 
of motion doesn't make contribution to the above distri- 
bution of constant-energy-time intervals. Second type of 
motion corresponds to the above mentioned synchronized 
motion. It is easy to show that the distribution for the 
constant energy time intervals for this case behaves as 
N(t) ~ t~ 3 . 

To explain this we note that in the synchronized regime 
billiard particle returns to its initial position during the 
kicking period. If in each initial condition we will shift its 
initial position to a infinitesimally small distance SL (see 
Fig. SJ its final (after one kicking period) position shifts 
to the same distance SL becoming closer to the kicking 
area. Then time t after which particle reaches the kicking 
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FIG. 4. Schematic picture 
area can be found from the following relation: 



where L is the distance between the particle's initial po- 
sition and the kicking area. 

Thus the number of initial conditions is related to the 
quantity SL as 

N(SL) = / p{q)d(q) (5) 
Jo 

where p is the density of particles. 

Since particles are uniformly distributed, p depends on 
q as 

p(q) = ii 2 (6) 

where 7 is coefficient of proportionality. 
Therefore it follows from the eq. ([5]) that 

N(SL)= 11 (SL) 3 (7) 

or, combining with the eq.@ we get finally: 

N(t) ~ N(SL) = 72 i- 3 (8) 

When considering not only a single trajectory, but 
rather an ensemble of particles, a diffusive growth of the 
ensemble averaged energy (Eb(t)) (averaged of 1000 tra- 
jectories) can be observed, see Fig. [SI More precisely, the 
energy growth linearly (normal diffusion) with the time 
t. The proportionality constant is naturally just the dif- 
fusion coefficient Db(a, T), so that (Eb(t)) — Db(a, T) ■ t. 
As indicated, Db depends on the coupling constant a and 
the kicking period T. Db increases monotonically with 
increasing a and decreases monotonically with increasing 
T, see Fig. 

It is reasonable to compare the diffusive growth of the 
ensemble averaged energy (Eb(t)) — Db(a,T) ■ t of the 



FIG. 5. Comparison of the time-dependence of the ensemble 
averaged energy of the kicked billiard and kicked rotor for 
a = 5.0, T = 1.0. 



FIG. 6. (Color online) The dependence of the diffusion coef- 
ficient Db on a and T for the kicked billiard. 



kicked billiard with the evolution of energy (E r (t)) of the 
kicked rotor (again with a being the coupling constant 
and T being the period between two successive kicks). 
From Fig[5j it can be seen that (E r {t)) — D r (a, T)-t, but 
with D r (a,T) <C D b (a,T), so the energy growths much 
faster in the case of the kicked billiard. We note that in 
the case of the kicked rotor the dynamics is effectively 
governed by a single parameter K r — aT only, see e.g 
[2j. The energy of the kicked rotor as a function of the 
dimensionless time n = t/T can than for K r > 5 be 
written as E r (n) — K^/A ■ n, so 

D r {a,T) = D r {K r ) = K 2 r /A. (9) 

In the kicked billiard the situation is different, the dy- 
namics depends on a and T individually, thus there is 
no simple representation for D b (a,T) as in Eq. (j^J pos- 
sible. The dependence of D b (a,T) on a and T is shown 
in Fig. [SI whereas D r (a,T) is shown in Fig. [BJ The 
large deviation between D b and D r becomes immediately 
clear when considering the maximum momentum trans- 
fer AP max at a single kick. In the case of the kicked 
billiard, AP max — S / na/a 1 whereas for the kicked rotor, 

AP max OL. 

The corresponding distributions iV r (AP) and N b (AP) 
(1000 initial conditions iterated until t — 10 6 T) are 
shown in Figs. 13 |9] (the subscripts r and b denote the 
kicked rotor, the kicked billiard respectively). Clearly, 
the available range in AP increases with increase of a 
in the case of the kicked billiard, and for a fixed a this 
range is much larger than in the case of the kicked ro- 
tor. Unlike the N r (AP), the distribution N b (AP) has a 
jump at AP = 0, this can be explained as following. The 
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FIG. 7. (Color online) The dependence of the diffusion coef- 
ficient D r on a and T for the kicked rotor. 

momentum transfer in the case of kicked rotor can be 
found from 2-D standard map as AP r = K r sm(0). For 
the kicked billiard the momentum transfer can be written 
as: 

AP = \P> -P\ = K b sin(6p) (10) 

where P and P', which are the momenta before and after 
the kick, respectively and b — 87r/a, K b = ba. 

Hence, for the momentum transfer distribution we get 

N b (AP) - N(p) ( arcsin( AP + ^ SP ) - arcsin(— ) ] 
V K K J 

(11) 

In the case of kicked rotor we have N r (9) — const, for 
K r > 5 (chaotic regime), while for kicked billiard as it 
was mentioned above, we have equally distributed kick 
occurrence over x, y. Therefore, the dependence of this 
distribution on p is given by Nip) — const\ + const2p. 
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It follows from the eq.(Q]| that AP = for p = and 
for p = a/4. However, the kick numbers at these two 
points are not equal, number of kicks at p = a/4 points 
is maximal, while for p — it is minimal that explains 
the difference between the shape of jV(AP) in Figs. [5] 
and El 




aP 



FIG. 8. (Color online) Distributions of momentum transfers 
AP for the kicked billiard for a = 5.0 and for T = 1.0; 
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FIG. 9. (Color online) Distributions of momentum transfers 
AP for the kicked rotor for a = 5.0 , T = 1.0. 

In Fig. I10[ the time-dependence of the average energy 
(E(t)) is plotted for different localizations of the kicking 
source and compared to the one when kicking source is 
located at the center of the billiard. The shift of the 
kick source leads to minor changes in (E(t)) only, so the 
dynamics is rather robust against this shift. We tested 
this for other delocalizations as well, with similar results, 
thus there are not shown here. 

So far, we have considered the kicked billiard with a 
potential which is very similar to the one of the kicked ro- 
tor, in particular with the same minimum and maximum 
values of the kicking potential. Despite this, however, 
the acceleration in the kicked billiard is much more pro- 
nounced compared to the kicked rotor, which is of course 
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FIG. 10. Time-dependence of the ensemble averaged energy 
for different localizations of the center of the kicking source 
in the billiard (a — 5, T — 1): The central localization is 
compared to the case, when kicking source is shifted to the 
left (position of the center at (0.6; 0.5)) and when it is shifted 
along the diagonal of billiard (center at (0.6; 0.6)). 

due to the fact that much larger AP are possible upon 
single kicks in the billiard. To make the two systems bet- 
ter comparable in terms of the momentum transfer, we 
modify the potential given in Eq. ([T]) slightly and obtain: 
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FIG. 11. Comparison of the time-dependence of the ensemble 
averaged energy of the kicked billiard and the kicked rotor 
with a = 5.0, T = 1.0. 

, ff^acos(^) V 5(t-nT) p<R 
V(x,y,t) = I 87r V a V (12) 

[0 p > R 

This potential leads to the same maximum and minimum 
momentum transfer as in the kicked rotor. Now, the 
energy of the kicked rotor growth faster than the one 
the kicked billiard, see Figs. [TTJ These reasons are the 
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following: 




of billiard walls, while in our case the energy grows due to 
the interaction with the kicking source. The latter mech- 
anism is close to the kicked rotor rather than to the bil- 
liard with time- dependent boundaries. Detailed study of 
Fermi acceleration in billiard geometries showed that the 
acceleration is possible for non-integrable billiards only 
However, the above studied kicked square bil- 



IAS 



Hard has integrable boundaries. Despite this acceleration 
is possible in this system. This is the main difference be- 
tween the above kicked billiard and other time- dependent 
billiards (with time-changing boundaries). 



FIG. 12. (Color online) The dependence of the diffusion co- 
efficient Dt on a and T for the kicked billiard. 

Firstly, in the kicked rotor at every t — nT, n = 
1,2,3,... the particle experiences a kick, a certain 
amount of momentum is transferred. In the kicked bil- 
liard this is not the particle gets a kick only if is at 
t = nT, n = 1,2,3,... inside the kicking area, see Fig. [T] 
Consequently, particles will experience fewer kicks during 
a certain time interval compared to the kicked rotor. 

We note that the above considered kicked billiard to 
some extent is equivalent to time- dependent Sinai scat- 
tering billiard. However, unlike usual Sinai billiard the 
scatterer is not static and switches on periodically. In ad- 
dition, in kicked billiard collision with the scatterer is in- 
elastic, while in Sinai billiard collision with the scatterer 
is elastic. Here we mention that time- dependent Sinai 
billiard has been studied earlier in the context of Fermi 
acceleration 
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However, these studies deal with 
the billiards whose boundaries are time-dependent, while 
in our system the boundaries of the billiard are static. 

Therefore exploring of such time- dependent billiards 
requires solving the classical equations of motion with 
time- dependent boundary conditions. Time- dependence 
in our system is caused by the kicking source and we do 
not need to consider time- dependent boundary conditions 
in the equations of motion. 

Also, as it was shown, under certain conditions particle 
motion in the kicked billiard can be localized (trapping of 
particle into the kicking source) in the kicking area. Such 
a trapping is not possible for usual time- dependent Sinai 
billiard. 

Finally, it should be noted that particle acceleration in 
kicked billiard can be considered as a kind of Fermi ac- 
celeration. In usual Fermi acceleration the growth of the 
kinetic energy of a billiard particle is caused by the motion 



le+06 




10000 



FIG. 13. (Color online) Time-dependence of the number of 
survived particles N(t) for the value a — 1.0 and calculated 
using eq. JT3). a = 0.001 x a, No = 10 6 ,po = 1.0. The inset: 
Distribution of the escaped particle momenta. 
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FIG. 14. (Color online) Time-dependence of the number of 
survived particles N(t) for different values of the delta-kick 
strength a.a = 0.001 x a, No = 10 6 ,p = 1.0. 

In this section we study particle transport in a kicked 
open billiard. The first system we will consider is a 
square billiard having one small hole with the size a. The 
type and location of the kicking source is the same as in 
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the previous section. Simulation of the system is per- 
formed for N non-interacting particles with randomly 
distributed initial positions, (x ,y Q ), anc ^ directions of 
the initial velocities, 0. The initial momenta of all par- 
ticles are assumed to be equal pq. Classical dynamics of 
unperturbed open billiards has been extensively studied 
both for integrable and chaotic geometries 24 1- [ill. in 
particular, it was found in the Refs. 24|- 26J that for the 
case of ergodic particle motion in an open billiard the 
particle's escape rate decreases exponentially, while for a 
non-chaotic system it decreases according to power law. 

In this work we explore behavior of the escape rate 
for the kicked open billiard. The quantity we want to 
calculate is the number of particles in billiard at the time 
t, N(t). The escape rate is related to N(t) as 

p{t) = 1 - N(t)/N , 

where No is the number of particles at t = 0. As it 
was shown recently 24, [2^], for chaotic billiards time- 
dependence of survival probability can be written as 



N(t) 

An 



exp(-Ci), 



(13) 



attached one incoming and two outgoing leads Fig. [15] 
Particles are assumed to come into the billiard from in- 
coming lead, while their escape from the billiard is pos- 
sible from all three leads. Reflection coefficient R is 
calculated as the escape rate from the incoming lead, 
while transmission coefficients T\ and T2 are the escape 
rates from the outgoing leads. In Fig. [16] these coeffi- 
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FIG. 15. (Color online) Kicked square billiard with three holes 



where 



c = 



with p being the absolute value of the billiard particle's 
momentum, S is the area of the billiard. Intuitively, 
one may expect similar decay in the case of kicked bil- 
liard (which is chaotic for arbitrary geometry of billiard 
boundaries). However, in this case momenta of billiard 
particles are not equal and p in eq. ([T3]) can be replaced 
by p' , most probable value of the escape momentum. In 
Fig. [32] we compare N(t) computed: 

i) using eq. ([T3]) and 

ii) by numerical modeling for three values of the kicking 
strength, a = 1.0. 

As it can be seen from these plots, increasing of the kick- 
ing force leads to decreasing of N(t) which is equivalent 
to increasing of the number of escaped particles. Such 
behavior can be explained by the fact that according to 
the previous section, for higher values of the kicking force 
the energy of the billiard particles growth linearly in time. 
Also, the plots in Fig. [14] show that the results of sim- 
ulation are in good agreement with the results obtained 
using eq. ([T3]) . 

To study transport properties of the kicked billiard sys- 
tem we consider square billiard having three holes with 



cients are plotted as a function of kicking strength. As 
it can be seen from these plots, the reflection coefficient 
is much smaller compared to transmission coefficients for 
the leads 1 and 2. This implies that such system can be 
used as a conductance amplifier. All the coefficients be- 
come a— independent for the values of a larger than 0.5. 
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FIG. 16. Transmission (T) and reflection (R) coefficients ver- 
sus kicking force a. a = 0.001 x a, No — 10 6 ,po = 1.0 
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IV. SUMMARY 



In this work we studied classical dynamics of a kicked 
particle whose motion is confined in a square billiard. 
The kick potential is considered as localized inside the 
billiard with central symmetric spatial distribution. It is 
found that for this type of kick potential the average en- 
ergy of the particle as a function of time grows faster than 
that of kicked free particle. This implies that the above 
explored system is more attractive from the view point 
of acceleration. Unlike the kicked free particle (kicked 
rotor) dynamics of the kicked billiard particle depends 
of perturbation parameters a and T separately. It also 
depends on the geometry of billiard and localization of 
the kick source. In addition, the case of the kick poten- 
tial that gives the same maximum and minimum values 
of momentum transfer as those for kicked rotor is con- 
sidered. It is found that in this system the acceleration 
is much weaker compared to the first system and even 
compared to kicked rotor. 

Also, we have studied particle transport in kicked open 
billiard with one and three holes. In particular, we cal- 
culated escape rate, transmission and reflection coeffi- 
cients. The latter quantities are calculated for a billiard 
with one incoming and two outgoing holes. Exponential 
decay (i.e., decreasing the number of billiard particle ac- 
cording to exponential law) of the system is found for the 
case of one-hole open billiard. It is shown for such bil- 
liard that increasing of the kicking force leads to rapid 
decay of the system. Transmission and reflection proba- 
bilities are studied for the case of a kicked billiard with 
three wholes. For this system injection of particles into 
the billiard occurs from one hole, while escape of par- 



ticles occurs from all three holes. It is shown that for 
such system transmission and reflection coefficients de- 
pend on the kick strength, a until we increase a up to 
certain value. Further increasing of a does not lead to 
changes in T and R. 

Finally, we note that the above considered kicked bil- 
liard is a kind of time- dependent billiards whose dynamics 
are completely different than that of static billiards. Re- 
cently billiards with time- dependent boundaries have been 
extensively studied }1 1] -J23l in the context of Fermi ac- 
celeration. In such billiards the dynamics is governed 
by the time-changing law of the boundaries. This time- 
dependence requires solving of classical equations of mo- 
tion with time- dependent boundary conditions. Unlike 
these billiards, the kicked billiard is perturbed by the kick- 
ing source located inside the billiard and the walls of the 
billiard are static. This fact causes the main difference 
in the dynamics of the kicked billiards and the billiards 
with moving walls. 

The importance of the above study is caused by its 
perspective relevance to the particle transport in various 
mesoscopic systems (e.g. as quantum dots, ratchets, nan- 
otubes etc.). Time-dependent external field can be used 
in these systems as an additional tool for manipulating by 
particle dynamics and transport. In this context an im- 
portant issue should be the extension of the above study 
for the case of corresponding quantum system. Such a 
study is currently in progress. 
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